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Abstract —This paper considers the problem of tracking a 
dynamic sparse channel in a broadband wireless communica¬ 
tion system. A probabilistic signal model is firstly proposed to 
describe the special features of temporal correlations of dynamic 
sparse channels: path delays change slowly over time, while 
path gains evolve faster. Based on such temporal correlations, 
we then propose the differential orthogonal matching pursuit 
(D-OMP) algorithm to track a dynamic sparse channel in a 
sequential way by updating the small channel variation over time. 
Compared with other channel tracking algorithms, simulation 
results demonstrate that the proposed D-OMP algorithm can 
track dynamic sparse channels faster with improved accuracy. 

I. Introduction 

In a broadband wireless communication system, channel 
state information (CSI) is required at the receiver due to the 
fact that the multipath fading channel distorts the received 
signals, especially when the channel is dynamically changing. 
Hence, accurate channel estimation and tracking become an 
important problem for communication over a dynamic wireless 
channel (lj. 

Various linear channel estimation methods with low com¬ 
putational complexity have been proposed 0, but their per¬ 
formance is often not robust enough to meet the requirement 
of communication systems with high rate and high mobility. 
Recently, a lot of physical channel measurements have verified 
that wireless channels exhibit sparsity, i.e., the dimension of 
a wireless channel may be large, but the number of active 
taps with significant power is usually small, especially in a 
broadband wireless communication system 0. By exploiting 
this channel sparsity, many nonlinear channel estimation meth¬ 
ods based on classical compressive sensing (CS) algorithms 
have been proposed to improve the estimation performance, 
such as orthogonal matching pursuit (OMP), compressive 
sampling matching pursuit (CoSaMP), and subspace pursuit 
(SP) 0-0. Compared with the conventional linear methods, 
CS-based channel estimation methods are able to achieve 
improved accuracy with reduced training resources 0. Further 
studies have uncovered additional channel characteristics, e.g., 
the temporal correlations of practical wireless channels: path 
delays have been shown to change slowly over time, while path 
gains evolve faster 0. However, the CS-based channel esti¬ 
mation methods ignore these temporal correlations of dynamic 
sparse channels and have to estimate them independently. 
By taking the temporal correlation into account, the adaptive 


simultaneous OMP (A-SOMP) algorithm has been proposed 
in 0 to obtain simultaneous multi-channel estimates based on 
the assumption that the dynamic channel estimates in several 
consecutive time slots share the same path delay set. However, 
the path delays of a dynamic sparse channel will change over 
time, or even there maybe some mutations, so the assumption 
in 0 is not always true in practice. Another attractive solution 
for tracking a dynamic sparse channel is the hierarchical 
Bayesian Kalman (HB-Kalman) filter based on Bayesian CS 
(BCS) ca, whereby the iterative re-estimation of the posterior 
covariance ECO is used to achieve accurate channel estimation 
when sudden changes happen to the dynamic sparse channel, 
but it suffers from slow tracking speed and high computational 
complexity. 

In this paper, we propose a dynamic CS algorithm called 
differential orthogonal matching pursuit (D-OMP) based on 
the standard OMP algorithm to track a dynamic sparse channel 
with fast tracking speed and low computational complexity. 
By exploiting the temporal correlation of a dynamic sparse 
channel, the proposed D-OMP algorithm only needs to de¬ 
tect the small variation of the dynamic sparse channel in a 
sequential way. Furthermore, an adaptive threshold based on 
the statistical analysis of the “equivalent” noise is proposed to 
accurately distinguish true, i.e., non-zero channel taps with low 
power from thermal noise. This is essentially different from the 
standard OMP algorithm in which the incorrect channel taps 
chosen in one iteration will never be removed, which finally 
leads to performance degradation. The performance analysis 
indicates that a non-zero channel tap can be detected with 
a high probability, while the noise in the estimate can be 
removed almost completely. Numerical simulations show that 
the proposed D-OMP algorithm can track dynamic channels 
faster and achieve more accurate channel estimates than other 
tracking algorithms. 

The remainder of this paper is organized as follows. The 
system model of dynamic sparse channels is described in Sec¬ 
tion II. Section III addresses the proposed D-OMP algorithm, 
together with the threshold based on noise statistics. Section 
IV presents the performance analysis, and simulation results 
are provided in Section V. Finally, conclusions are drawn in 
Section VI. 

Notation: We use upper-case and lower-case boldface letters 
to denote matrices and vectors, respectively; (-) T , {-) H , (-) -1 , 



Fig. 1. Illustration of the dynamic Vehicular B channel with a velocity of 
120 km/h in three consecutive time slots. 

(•)', and || • ||p denote the transpose, conjugate transpose, 
matrix inversion, Moore-Penrose matrix inversion, and l p norm 
operation, respectively; |Ij denotes the number of elements in 
set 1’; hp denotes the entries of the vector h in the set 1'; <t>p 
denotes the submatrix comprising the T columns of <f>. 

II. System Model 

In this paper, we consider the problem of tracking a dynamic 
sparse channel {h^}^ from under-sampled noisy measure¬ 
ments {y^}tLi in T time slots, in a broadband wireless com¬ 
munication system. The single time slot measurement vector 
y(*) = [y^\y • ,y^)] T is usually obtained through the 
linear measurement process 

y(‘) = $(f) h (‘)+n (<) , t = 1, 2, ■ ■ ■ , T, (1) 

where hW = [h^\ h^\ ■ ■ • ,denotes the time-domain 
discrete channel vector with N taps, where N > M, and 
the additive white Gaussian noise (AWGN) is modeled as 
~ CJ\f(0, where I m denotes the identity matrix 

of size M x M. The measurement matrix of size M x N 
is a Toeplitz matrix determined by the training sequence 
c = [co, Ci, • ■ ■ , cm] T in the time or frequency domain, e.g., 
the well known pseudo noise (PN) sequence used in CDMA 
systems 0. Thus the measurement matrix <!>(*' becomes 
time-invariant 0, he., = 3?^ = ••• = = 

# = [4> 1 , <f> 2 , • ■ • ,4>n]- Since the training sequence is usually 
designed in a random way, i.e., c* ~ CJ\f(0, a ;?), we can obtain 
that </>fj ~ CAf(0, cr^), where cr^ = and <denotes 
the element of the matrix <j>W. Furthermore, based on this 
Gaussian assumption, the columns of the $ matrix are semi- 
orthogonal, i.e., 4>Y<pj S3 0,i Y 3- 

The temporal correlations of practical wireless channels 
have been verified through analysis and experiments, even 
when the channels are varying fast 0- Fig. m illustrates the 
time-domain impulse response of the dynamic Rayleigh fading 
Vehicular B channel with a velocity of 120 km/h in three 


consecutive time slots EH. It is clear that path delays of 
such dynamic sparse channels change slowly, while path gains 
evolve faster. In order to characterize the temporal correlations 
of a dynamic sparse channel, we adopt a probabilistic signal 
model with two time series vectors {sW}^ and {aS^}JLi, 
where the binary vector • • • , s^] T is used 

to describe the temporal evolution of the channel paths, 
while the complex-valued vector = [a/" 1 , o.if 1 , ■ • • , a^] T 
characterizes the temporal variation of the path gains. Hence, 
the dynamic sparse channel can be modelled as 

=4 t)a ?\ t = 1,2,--- ,T, l<i< N (2) 

where the binary sf' 1 € {0,1} denotes whether there is a non¬ 
zero channel tap at index i in the / th time slot, and a/ * denotes 
the corresponding path gain. Thus, the path delay set of 
the dynamic sparse channel can be represented as = {i : 
= !}• 

We model the time series {s^}Y—i as a discrete Markov 
process with two transition probabilities E): Po->i = 
P{s[ t] = l|sf _1) = 0} and pi^ 0 = P{s = 0|s- 4-1 ^ = 1}. 
Without loss of generality, {s| 1 '*}A =1 is initialized as indepen¬ 
dent Bernoulli random variables based on the probability pi, 
i.e., ^Bernouli(pi), i = 1, 2, • • ■ , N. Similarly, we model 
the time series {a^}f =1 as a Gauss-Markov process, and each 
path gain evolves independently as 

af ] =af _1) +wf ) , t = 1,2,--- ,T, l<i< N (3) 

where wW = w^\ ■ • • ,w^] T ~ CAf(0, con¬ 

trols the temporal correlation of the path gains G3). We model 
the initial distribution of path gains as a Gaussian distribution, 
i-e., ~C/V(0,f7^Iiv). 

The probabilities po->i an d pi-j-o are usually small to model 
the slow changing of the path delays, while {a^}tLi changes 
in each time slot to model the faster changing of the path gains. 
Furthermore, the aforementioned parameters will be given in 
detail in Section V. 

III. Proposed D-OMP Algorithm 

In this section, the proposed D-OMP algorithm is explained 
in detail. Then the threshold used in the algorithm is derived 
based on noise statistics. 

A. D-OMP Algorithm 

We propose the D-OMP algorithm to track a dynamic sparse 
channel, which obtains the final estimates in a sequential way 
by updating the small variation of the dynamic sparse channel. 
The key idea of the proposed D-OMP algorithm is that, the 
major information of the channel in the current time slot can 
be obtained from the estimation results in previous time slots 
due to the temporal correlations of the dynamic sparse channel, 
and then the small variation of the channel can be estimated 
with low complexity to refine the final estimate result. This is 
essentially different from the standard OMP algorithm which 
ignores any temporal correlation of the channel taps 0. 
The proposed algorithm also differs from the HB-Kalman 











algorithm utilizing the iterative re-estimation of the posterior 
covariance to obtain channel estimation, which leads to a slow 
tracking speed and high computational complexity. 


Algorithm 1 D-OMP Algorithm 

Input: 

Received signals: {y^}^; 

Measurement matrix: $ = [</> 1; <f> 2 , ■ ■ • , 

Threshold: P t h- 

Output: 

Channel estimates: {iv 4 )}^. 

1: Initialization : 

2: AW = 0, h(°) = 0. 

3: for t = 1 to T do 

4: y+ = y(‘) - <Ph (4-1) . 

5: A+ = argmax^A, = \4>f y + | : i £ A^" * 1 )}. 

6: A (4) = A (4_1) U A+. 

7: = argminhlU^h — y ^ 4 ^||2 : supp(h) C A^ 4 )}. 

8: A~ = {i:\hf ) \<P th ,ieA^}. 

9: = 0, A^ 4 ) = AW - A - . 

10 : end for 

11: return {hW}^. 


The pseudocode of the proposed D-OMP algorithm is 
provided in Algorithm 1, and each loop is comprised of the 
following three parts: 

1) Incremental Support Detection (step 4^6): Firstly, 
we obtain the incremental observations y + by subtracting 
<bh( 4-1 ) from yW in step 4, where If 4 1 ) denotes the estimate 
of lA 4-1 ). Here we use 4>lr 4 ^ instead of y( 4 ^ to avoid 
the impact of the noise if 4-1 A Different from the way used 
in the standard OMP algorithm which iteratively detects the 
candidate supports by selecting the element that correlates 
best with the residual signal 0, here we only detect the 
incremental support A + in step 5, which is able to acquire 
the appearing channel tap as sf -1 ^ = 0 —> s = 1. Then, 
the temporary estimate of the path delays set A W j s obtained 
by merging A ft_ l * which contains the persistent channel taps, 
i.e., i £ A^ 4-1 ) D AW, and A + which contains the appearing 
channel tap. 

2) Channel Estimate Update (step 7): Step 7 aims to update 
the channel estimate based on the temporary estimate of the 
path delays set AW. Intuitively, the variation of the channel 
estimate is comprised of three parts: a) The appearing channel 
tap af* when sf -1 '* = 0 —> s'f 1 = 1; b) The disappearing 
channel tap —d* ^ when s^ 1 ’ ) = 1 —> s^ = 0; c) 
The smooth variation of the path gains wf\i £ AW). The 
minimization problem in step 7 corresponds to the main 
computational burden of the D-OMP algorithm, which can 
be realized by the standard least-squares (LS) technique with 
low complexity m. 

3) Disappearing Channel Taps Removal (step 8~9): When 

a tap disappears as sf 1 3 ’ > = 1 —> s,W = 0, the 

estimated tap in the previous (t — l)th time slot will still remain 
in the estimated path delay set of the current time slot since 


AW = A 1 - 4 ' WjA^. The path gain of this disappearing channel 
tap A-*" 1 will be small. Therefore, we propose a threshold 
P t h to judge whether the nonzero elements in hW are non¬ 
zero channel taps or noise in step 8. After removing the 
disappearing channel taps, the final channel estimate can be 
obtained in step 9. 

Unlike the standard OMP algorithm which needs to solve 
the minimization problem for many times according to the 
sparsity level of the wireless channel 0, the proposed D- 
OMP algorithm only needs to solve the minimization problem 
once, so the computational complexity can be dramatically 
reduced. Furthermore, unlike the hard threshold 0 or the 
rough criterion 0, we rely on the noise analysis to select 
the threshold P t h, which will be derived in the following 
subsection. 

B. Threshold Based on Noise Statistics 

Here we consider a single time slot (we omit the time slot 
superscript t), and rewrite © as 

y = <bh + n = 3>(h + <b^n) = $(h + n ), (4) 

where denotes the pseudoinverse of 3?. 

Generally, linear channel estimation methods are designed 
to directly approach h + n , so their performance is limited 
by the contamination of the “equivalent” noise n = ' n in 

(4). The advantage of CS-based nonlinear estimation methods 
is that they reconstruct the sparse channel h and avoid part of 
the interference caused by the “equivalent” noise n . Since the 
basic idea of the greedy CS methods is to iteratively search 
the channel taps from the strongest one to the weakest one 
0, 0, it is difficult for these algorithms to judge whether 
the searched taps are true, i.e., non-zero channel taps or noise, 
when they are weak. 

In order to further reduce the interference caused by n 
under the CS framework, we design the threshold P t h based 
on the statistical properties of the noise n for the pro¬ 
posed D-OMP algorithm to separate non-zero channel taps 
from noise as much as possible. According to ©, we have 
n = The mean of the vector n can be easily obtained 
as E{n } = E{$^}E{n} = 0 due to the fact that the 
measurement matrix $ and the noise vector n are independent. 
Then, the covariance of the elements in equivalent noise n 
can be derived as 

Var{rij} = Var{0|n} = MNar{(j)\ -}\as:{nj} 

= MolN ar{<4}, (5) 

where (j)\ denotes the Ah row of the matrix <!>' and cp[ is the 

t ^ 

yth element of <p). Furthermore, we have 

= ~ ( 6 ) 

where we use No^Im to approximate <!><!> 11 due to the 
randomness of the matrix 4>. Then we can obtain the variance 







of rij as 


Var{n'} 


{Naif 


Var« 7 } 


tvv|- 


(7) 


Finally, the statistical properties of the elements in n can be 
represented as 


Hi ~ CAf{0,a 2 ,) 


Cj\f (o 

1 ’ jwr - 


( 8 ) 


With the help of the statistical properties of the noise vector 
n , the threshold P t h can be selected as 


Pth = OLO n ' 


a.o n \[M 

a^N 


(9) 


where the variance a 2 of the noise vector n can be usually 
obtained at the receiver na. The coefficient a can be set as 
a = 3 so that the strength of the noise element | n t | will be 
smaller than Pth with high probability of 99.73%, based on 
the assumed Gaussian distribution. 


IV. Performance Analysis 


A. Persistent Channel Tap Detection 

In section III-B, the proposed threshold Pth = otcr n ' is 
used to remove the disappearing channel taps from the channel 
estimate result with high probability. On the other hand, as the 
proposed scheme may also remove a persistent channel tap by 
mistake, the correct detection probability of the persistent tap 
for the scheme is also important to guarantee its performance 
in practice, which will be analyzed below. 

Firstly, the SNR 7 at the receiver can be obtained ma and 
represented as 


7 = lOlog 


Var{*I>h} 

Var{n} 


lOlog 




lOlog 


MKai 

N 2 a 2 , 


, ( 10 ) 


where K = Npi denotes the sparsity level of the dynamic 
sparse channel. Thus, the threshold P t h can be obtained as 


Pth = aa n ' 


aoh 


MK 

IQtoN 2 ' 


( 11 ) 


After the threshold has been obtained, we can calculate the 
detection probability of a persistent channel tap as P(h [G 
A), which can be derived by using the normal Gaussian 
distribution probability function T as 


p(h™ e A) 


P{\hf\>p th ) 

2 ( 1 - tf(a 


■A AT?. > 


IOtb N 2 ' 


( 12 ) 


For example, in a typical wireless communication system with 
M = 200, N = 400, K = 5, and 7 = 15 dB, the probability 
of a persistent tap being larger than P t h is about 96.64%. 


B. Appearing Channel Tap Detection 

As the D-OMP algorithm is proposed to track a dynamic 
channel, the condition of appearing channel tap detection is 
important to be analyzed. For the sake of a simplified analysis, 
we assume that the channel estimate in time slot (t — 1 ) 
is accurate without loss of generality, i.e., 

When there is an appearing channel tap from time slot (f — 1) 
to (f) (for simplicity, we assume other channel taps stay 
unchanged, i.e., hj^ = hj,* 1 \vfc € A^ _ 1 l), the step 4 in 
Algorithm 1 can be rewritten as 

y+ = yw _ $ h (*-l) = 0,/jW + n, (13) 

where hf' 1 (i ^ A*7 -1 )) denotes the appearing channel tap at 
the time index i during the /th time slot of the dynamic sparse 
channel. 

According to the detection rule in step 5 of Algorithm 1, 
the appearing channel tap lip is successfully detected if and 
only if 

A j < A i, Vj ^i,j i A (t_1) , (14) 


where A, can be derived as 


A* = \4>fy + \ = l0fW>A w + n )l > 

Then Xj can be similarly derived as 


I0fn|. 

(15) 


A j = \<pf y + | = 1 4>f + n)| < 1 4>f (friWh^l + 1 4>f n|. 

(16) 

Using ®. we can obtain the condition under which the 
appearing channel tap is successfully detected: 


1^1 > 


|0f n| + |0f n| 

\<t>?<t>i\ - \4>f 4>i\’ 


(17) 


which means the path gain of the appearing channel tap should 
be larger than a minimum value. 

Intuitively, the idea behind the incremental support detection 
of the proposed D-OMP algorithm is to find the elements A + 
that correlate best with the incremental observations y + , but 
the previously detected support A'* -1 ' will not be considered. 
Due to the fact that the appearing channel tap can be detected 
immediately if it is larger than a certain value derived in 
(fT71l . while the posterior covariance matrix of the HB-Kalman 
method ifTOil needs several time slots to converge, the condition 
of incremental support detection of the proposed D-OMP 
algorithm is much easier to satisfy than that of the HB-Kalman 
method. 


C. Computational Complexity 

The computational complexity of the proposed D-OMP 
algorithm in terms of the required number of complex multi¬ 
plications includes the following three parts: 

1) In the incremental detection part, the complexity is 
O(NM) for the calculation of $h( t_1 ) and \<j>fy + \. 

2) In the channel estimate update part, the LS problem 
argminh{||$h — y + ||2 : supp(h) C A^} can be 
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Fig. 2. MSE performance comparison against time slot t. 

implemented with the 0(MK 2 ) complexity by using 
the Gram-Schmidt algorithm m. 

3) In the disappearing taps removal part, the complexity 
is O(M) for the comparison of taps with the threshold 

Pth- 

To sum up, the total complexity of the proposed D- 
OMP algorithm is 0{M{N + K 2 + 1 )T) to track a dy¬ 
namic sparse channel over T time slots. Compared with 
the 0(KM(N + K 2 )T) complexity of the standard OMP 
algorithm, the complexity of the proposed D-OMP algorithm 
is reduced approximately by a factor of K. Considering 
the re-estimation of the posterior covariance matrix with J 
iterations in each time slot ED, the computational complexity 
of HB-Kalman filter (TO) is 0(JKM(N + K 2 )T), which 
is much higher than the complexity for both the standard 
OMP algorithm and the proposed D-OMP algorithm. The 
complexity reduction comes from the fact that the proposed 
D-OMP algorithm obtains the major information of the current 
channel from the estimation results in previous time slots, and 
then refines the final estimation result by detecting the small 
variation of the channel. 

V. Simulation Results and Discussion 

This section investigates the performance of the proposed 
D-OMP algorithm compared with the conventional linear 
algorithm 12, standard OMP algorithm ( 6 ), A-SOMP algo¬ 
rithm 0, and HB-Kalman algorithm uni. The parameters 
mentioned in Section II (system model) are set as follows: 

1) The size of the matrix $ is TV = 400 and M = 200; 

2) The probability pi is set as 0.025, which means the 

average channel sparsity level is K = Np± = 10; 3) A 
channel tap will appear or disappear over every l/(iCpi_>o) = 
l/((iV — K)po-yi) = 10 time slots on average corresponding 
to probabilities pi_> 0 = 0.01 and p 0 -n = PiPi-s-o/(l — Pi); 
4) a a = 0.05, er„ = 0.05, a h = 1, = 1. 

Fig. m shows the mean squared error (MSE) performance 
against time slot t for the five channel estimation methods 



Fig. 3. The correct detection probability comparison against SNR. 



Fig. 4. MSE performance comparison against SNR. 


mentioned above. It is clear that A-SOMP and HB-Kalmam 
achieve lower error levels than the standard OMP algorithm, 
while the conventional linear method performs worst. The 
MSE performance of the proposed D-OMP algorithm is the 
best, as the temporal correlations of the dynamic sparse 
channel are efficiently exploited. More importantly, when two 
channel taps suddenly disappear in time slot t = 80, the HB- 
Kalman algorithm requires about 20 time slots to detect this, 
while the proposed D-OMP algorithm is able to detect the 
change immediately, which confirms the tracking capability 
of the proposed scheme as discussed in section IV-B. 

Fig.[3]shows the correct detection probability of a persistent 
channel tap against SNR. It is evident that the hard threshold 
used in many CS-based channel estimation methods si is 
not adapted to the SNR. The rough criterion (9J can im¬ 
prove the correct detection probability by using the statistical 
information of the channel. For the proposed threshold P t h 
based on noise statistics, the correct detection probability can 
be improved further, which is close to the theoretical bound 
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derived in (ITll i. 

Fig. a shows the MSE performance comparison against 
SNR for the five channel estimation methods. It is clear that 
the standard OMP algorithm outperforms the linear method 
by about 1 dB, where the benefit comes from utilizing the 
channel sparsity. Further, A-SOMP and HB-Kalman are better 
than the standard OMP algorithm by about 2 dB, since they 
partially consider the temporal correlations of the dynamic 
sparse channel. For the proposed D-OMP algorithm, it is 
evident that another 2 dB SNR gain can be achieved due to 
its capability to track the dynamic sparse channel rapidly and 
detect the non-zero channel taps accurately as discussed in 
Section IV. 

VI. Conclusion 

In this paper, we have proposed a novel dynamic CS 
algorithm called D-OMP to rapidly track a dynamic sparse 
broadband communication channel in a sequential way by 
updating only the small channel variation over time. An 
adaptive threshold based on noise statistics is used in the 
proposed D-OMP algorithm to remove the disappearing chan¬ 
nel taps. Moreover, the performance analysis provides the 
correct detection probability of persistent channel taps as well 
as the condition of appearing channel tap detection for the 
dynamic sparse channel with low computational complexity. 
Finally, simulation results demonstrate that the proposed D- 
OMP algorithm is able to rapidly detect appearing channel 
taps, and achieves about 2 dB gain compared with the recently 
proposed dynamic channel tracking algorithms. 
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